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Abstract 


In this study, a fractal model is developed to predict the permeability and liquid water relative permeability of the GDL (TGP-H-120 carbon 
paper) in proton exchange membrane fuel cells (PEMFCs), based on the micrographs (by SEM, i.e. scanning electron microscope) of the TGP- 
H-120. Pore size distribution (PSD), maximum pore size, porosity, diameter of the carbon fiber, pore tortuosity, area dimension, hydrophilicity or 
hydrophobicity, the thickness of GDL and saturation are involved in this model. The model was validated by comparison between the predicted 
results and experimental data. The results indicate that the water relative permeability in the hydrophobicity case is much higher than in the 
hydrophilicity case. So, a hydrophobic carbon paper is preferred for efficient removal of liquid water from the cathode of PEMFCs. 


© 2006 Elsevier B.V. All rights reserved. 
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1. Introduction 


Carbon paper is commonly used as the gas diffusion layer 
(GDL) in the application of PEMFCs, which is composed of 
randomly packed carbon fibers. The GDL is an important part of 
PEMFC s. It has several functions: reactant permeability, product 
permeability (i.e. liquid water removal), electronic conductiv- 
ity, heat conductivity and mechanical strength for providing 
mechanical support to the MEA. The properties of the GDL 
directly influence the performance of PEMFCs [1-5], which 
are determined by the GDL’s microstructure. However, there 
are only a few reports on the relationship between the GDL’s 
properties and microstructure. Newman and co-workers [6] 
developed a statistical method for predicting the pore size dis- 
tribution (PSD), permeability and capillary pressure in a GDL. 
In the present work, a fractal model is developed for predict- 
ing the absolute permeability and water relative permeability 
based on the microstructure of the GDL obtained by SEM. 
Pore size distribution, maximum pore size, porosity, diameter of 


* Corresponding author. Tel.: +86 411 88993626. 
E-mail address: zczhao55 @ 163.com (Z. Zhao). 


0378-7753/$ — see front matter © 2006 Elsevier B.V. All rights reserved. 
doi:10.1016/j.jpowsour.2006.09.054 


the carbon fiber, pore tortuosity, area dimension, hydrophilicity 
or hydrophobicity, and the thickness of GDL as well as sat- 
uration are involved in this model, which provides an insight 
into the properties of the GDL and the water management in 
PEMFCs. 


2. Fractal theory 


Euclidean geometry describes ordered objects using inte- 
ger dimensions 0, 1, 2 and 3, respectively. However, numerous 
objects found in Nature [7] such as rough surfaces, coastlines, 
mountains, and islands, are disordered and irregular, and they do 
not follow the Euclidean description since their length, area and 
volume are scale-dependent [8]. Such objects are called frac- 
tals, and are described using a non-integral dimension called the 
fractal dimension [8]. The measure of a fractal structure, M(L) 
is related to the length scale, L, through a power law in the form 
of [7]: 


M(L) L”! (1) 


where ‘~’ means ‘scales as’, and Dr is the fractal dimension 
of the structure. Associated with Eq. (1), the property of self- 
similarity implies that the value of Df calculated from the rela- 


G. He et al. / Journal of Power Sources 163 (2007) 846-852 847 


Nomenclature 

Ao unit area (cm?) 

d fiber diameter (um) 
di tortuosity dimension 


Dg area dimension 

fa fraction of hydrophilic pores 
K permeability (cm?) 

ls length of space (cm) 

Lo straight length (cm) 

p pressure 

q flow rate in a capillary 

Q total flow rate 

s saturation 


Greek letters 


E porosity 

0 contact angle 

À pore diameter (cm) 
o surface tension 
Subscripts 

c critical 


cap of capillary 
max maximum 


min minimum 

TW of water relative 

rwi of water relative for hydrophilicity case 
rwo of water relative for hydrophobicity case 


tionship in Eq. (1) remains constant over a range of length scales, 
L. Geometric constructs such as the Sierpinsky gasket and the 
Koch curve are the examples of exact fractals which exhibit 
identical self-similarity over an infinite range of length scales. 
However, self-similarity in a global sense is seldom observed in 
nature and the fractal description is based on a statistical self- 
similarity, i.e. the objects exhibit self-similarity in some average 
sense, over a certain local range of length scale, L, relevant to the 
problem [9]. The concept of statistical self-similarity introduced 
above is used in the following section to develop a geometric 
description of the GDL in PEMFCs. 


2.1. Fractal characteristics of GDL 


Carbon paper is mainly used as the GDL of PEMFCs, whichis 
composed of randomly packed carbon fibers [10]. It has numer- 
ous pores with various sizes in the through-plane direction, and 
can be considered as a bundle of tortuous capillary tubes with 
variable radius for the two-dimensional case. Let the diameter 
of a capillary tube in GDL be å, and its tortuous length along 
the flow direction be L(A). The relationship between them also 
exhibits the fractal scaling law [11]: 


LO toy 
(a) ° 


where Lo is the representative length of a straight capillary, and 
dı is the tortuosity dimension, with 1 <d, <2. Large value of 
d within this range corresponds to a highly tortuous capillary, 
while d= 1 denotes a straight capillary pathway, d;=2, corre- 
sponds to a highly tortuous line that fills a plane [12]. 

Besides the convolutedness of the capillary pathways, the 
relationship between the number of pores and the pore size A is 
another important property. The pores in a porous medium are 
analogous to the islands or lakes on the earth [13]. The cumu- 
lative size distribution of them follows the power law relation 
[13]: 


À De 
N(L > 1) = (==) (3) 


where N(L > à) represents the total number of pores with diam- 
eter greater than A, Amax is the maximum pore diameter. From 
Eq. (3), the total pore number can be obtained 


Amax Dt 

N,(L = Amin) = ( ) (4) 
Amin 

and 

—dN = Dea Pt AT OPHD an (5) 


The negative sign in Eq. (5) implies that the pore population 
decreases with the increase of pore size, and dN>0. Eqs. (1) 
and (3) hold true for both exactly and statistically self-similar 
fractal geometries. Here, 1 < Df < 2, for two-dimensional space. 
According to Eqs. (4) and (5): 

-< = Dp ATOD dà = f(A) dà (6) 


t 


where f(A) is the probability density function, and 


CO Amax Am; Dt 
J fa)da = f fda =1— (=) (7) 
—0o Amin max 


The above equation equals unity if and only if 


Dees 
min _ 0 (8) 


Amax 


It implies that Amin K Amax must be satisfied for fractal analysis 
of a porous media, otherwise the porous media is a non-fractal 
medium. As for the GDL of PEMFCs, it holds true according to 
the measured pore size distribution [10], so it is a fractal. 


2.2. Fractal model for permeability of GDL 


The total volumetric flow rate, Q, through a unit cell is a sum 
of the flow rates through all the individual capillaries. The flow 
rate through a single tortuous capillary is given by modifying 
the well-known Hagen-—Poiseulle equation [14] to give 

m Ap ìf 


qa) = T28 LO) T (9) 


Then, the total volumetric flow rate in a unit cell is 


Amax 
Q= -f qa) dN(A) (10) 


min 
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The permeability can be expressed as follows according to the 
Dacy’s law: 


1— (Amin Amad i 


L 
po POS T na ip; 5 
ApAg 128 3 +d — Df 
T .3+d 7 -1-d 1 
= —_p3tay ld De 11 
pee "34d — Di 


where Ap = Li is the unit area. 

Thus, Eq. (11) gives the absolute permeability of carbon 
paper obtained by the fractal theory, and the maximum pore 
size (Amax), the thickness of the carbon paper (Lo =0.037 cm 
[15] for TGP-H-120), the tortuosity dimension (d+), and the area 
dimension (Df) are involved. 


2.3. Fractal model for water relative permeability of GDL 


Capillary-driven liquid water behavior in GDL is related 
to the wetting property of GDL for the two-phase flow situa- 
tion. First, for hydrophilicity case, water transports by wick- 
ing function, and the capillary force can be calculated as 
Pcap =40 cos(0)/Àc, so the wicking function increases while pore 
size decreases. It means that water prior to transports through 
the small pores and the critical pore diameter A, (i.e. the maxi- 
mum size of the pore that water can flow through) increases as 
saturation increases. The relationship between critical diameter 
and saturation is 


Àc 2 Df »—Dr-1 
i Jie, A? DADE ATPM! da 
T SAPS A2 DADE AP] dd 


Àmin 


(12) 


According to Eqs. (10) and (11), the water relative permeability 
can be determined as 


„ _— KLoQw/ApAo _ Qw _ Sisin (ANC) 
"i aLoQ/ApAo Q — [E*N 


Àmin 


(13) 


Submitting Eqs. (5), (8) and (9) into the equation above, then 


3+d,— Df 3+d— Df 3+d,— Df 
baa BAD BE ( 2) N 
wi — a 3d Df — 9 3+a—Dp 

Amar! j Amin — Amax 


Eq. (14) is the final expression for liquid water relative perme- 
ability of the GDL for the hydrophilicity case, where A, can be 
obtained in Eq. (12). 

As for the hydrophobicity case, initial pressure is needed 
for water transport due to the hydrophobicity, which can also 
be calculated as Pcap =40 cos(@)/Ac. It means the lager pores 
need smaller initial pressure. Water prefers moving through large 
pores, and the critical pore size (i.e. the minimum size of the 
pore that water can flood) decreases as saturation increases, so 
the saturation can be determined as 


fi A? DA Pt ATP! dà 


d 
= Fi 15 
T GREE AE Dea PE APT a Ha 
The water relative permeability is 
_ pLoQw/ApAo_ Qw _ Sa) ANO) 5 


kiwo = = = 
M pLoQ/ApAo Q — [q dN) 


Submitting Eqs. (5), (8) and (9) into equation above, then 


y3ta—Dp _ }3+di— Df À 3+d;— Df 
kwo = Fy eam, = 1- (=) g 
À x! t= Nen Amax 


ma: 


Eq. (17) gives the final expression for liquid water relative per- 
meability of GDL for the hydrophobicity case, where A, can be 
obtained in Eq. (15). 

However, the carbon paper is treated to be waterproof in the 
actual application to PEMFCs. There may be two kinds of pores 
in a GDL, i.e. hydrophilic pores and hydrophobic pores coex- 
isting. In this case, the water relative permeability is related to 
the fraction of the hydrophilic (or hydrophobic) pores to the 
total pores. far (or fuo) and the value of water relative perme- 
ability can be obtained by averaging the relative permeability 
of the hydrophilicity case and the relative permeability of the 
hydrophobicity case [6]: 


krw = Sutkrwi +- SaDkrwo (18) 
2.4, Determination of Df 


The area dimension Df can be determined by the box- 
counting method [16]. This method is based on the image anal- 
ysis of a unit cell or a sufficiently large cross-section of a sample 
along a plane normal to the flow direction. In this method, the 
cross-section under consideration is discretized using square 
boxes with size A, and the number of boxes, M(A) required to 
completely cover the pore areas is counted. The area dimension, 
Dp can be determined by the value of the slope of a linear fit on 
a logarithmic plot of the cumulative number of pores versus the 
pore size À. 

Fig. 1 shows a micrograph (by SEM) of the GDL (TGP-H- 
120), where the carbon fibers and the pores formed by random 
packed fibers can be easily distinguished. Fig. 2 shows the 
data obtained by box-counting method with the length scale 


Fig. 1. A micrograph of TGP-H-120 (through-plane direction). 
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log N(A) 


D=-slope=1.8991 


0.5 1.0 1.6 2.0 
log(A) 


Fig. 2. Schematic illustration of the procedure for evaluation of the tortuosity 
dimension, Df. 


range from 2 to 60 pixels. The data follows a closely linear 
relationship on the logarithmic scale, and this confirms the 
statistical fractal nature of the microstructure of GDL. The 
value of area dimension is 1.8991 according to the slope of the 
data line. 


2.5. Determination of tortuosity dimension, d; 


As mentioned previously, the tortuosity dimension represents 
the extent of convolution of the capillary pathways for water 
flowing through GDL. Since the tortuosity dimension pertains 
to the capillary pathways in the direction of the overall flow for 
the permeability evaluation, the microstructural images consid- 
ered for the calculation of d, are those of the GDL cross-section 
parallel to the flow direction (along the direction of thickness). 
Furthermore, since the tortuosity of the flow pathways results 
from the convolution of the boundaries of the porous regions 
as seen in a two-dimensional GDL cross-section, the tortuosity 
dimension may therefore be evaluated as the fractal dimension of 
the perimeter of the porous regions, which may also be obtained 
by the box-counting method mentioned above. Fig. 4 shows a 
log-log plot of M(L) against L obtained by the box-mounting 
method for the micrograph of TGP-H-120 shown in Fig. 3. As 
expected, the number of boxes required to cover the perime- 
ter decreases with increasing box size, and furthermore, the data 
points closely follow a straight line in the log(M(L))—log(Z) plot. 
The value of dt is 1.2507 according to the slope of the linear fit 
in Fig. 4. 


2.6. Determination of Xmax 


The maximum pore size appeared in the permeability expres- 
sion corresponding to the pore space formed between the carbon 
fibers in GDL. The spaces between carbon fibers are a func- 
tion of the architectural parameters of the GDL, and the GDL 
considered in this study consists of random packed ideally equi- 
spaced aligned carbon fiber screens [17] shown in Fig. 5. For 
this structure, the maximum pore space may be considered to 


Fig. 3. A micrograph of the GDL cross-section (TGP-H-120). 
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Fig. 4. Schematic illustration of the procedure for evaluation of the tortuosity 
dimension, d. 


Fig. 5. An ideal carbon fiber screen. 
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Fig. 6. A micrograph of the carbon fiber in GDL (TGP-H-120). 


be resulted from a single ordered carbon fiber screen shown in 
Fig. 5. According to this equation 


12 

c= s 
(ls +d}? 
Here, /, is the length of the space, and d is the diameter of a 
carbon fiber. According to the micrograph shown in Fig. 6, it 
possess the value of 12.5 um, and the porosity of TGP-H-120 is 


0.78 (Toray Carbon Paper properties [10]), then the maximum 
pore size can be 


(19) 


dena = 4 / L fae (20) 
3. Results and discussion 
3.1. Predicted permeability 


According to Eqs. (11), (19) and (20), the permeability of 
TGP-H-120 can be obtained. The permeability can also be 
obtained by the well-known Kozeny—Carman equation: 


ed? 


K = —___ 21 
16K,(1 — £)? 2D 


For fibrous media, Kx is taken as 6 [17]. The predicted result 
by Kozeny—Carman equation and the fractal model in this study 
are compared with the experimental data by Williams et al. [1] 
(TGP-H-120, the value of permeability is 8.69E—12 m7). The 
result of the present fractal model fits the experimental data well 
as shown in Fig. 7, so, the fractal model points out a new way 
to estimate the permeability of GDL in PEMFCs. 

Fig. 8 shows the permeabilities for different tortuosity dimen- 
sion, which reflect the extension of tortuosity in GDL. It can be 
seen that the permeability increases while dt decreases. When 
d,= 1, the permeability have the maximum value, which is con- 
sistent with the theory mentioned above, i.e. while d,=1, the 
capillaries are straight. It is advantageous for water flowing 
through. Fig. 9 shows the permeabilities for different area dimen- 
sion Df. It can be seen that the permeability increases with 
increase of Df. 
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Fig. 7. A comparison between the fractal model predicting permeability and test 
data. 
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Fig. 8. Permeability as function of d. 


3.2. Predicted water relative permeability 

Fig. 10 shows the water relative permeabilities for the 
hydrophilicity case, hydrophobicity case, fg = 0.3, and 0.8. The 
results indicate that water relative permeability of the hydropho- 
bicity case is much lager than that of the hydrophilicity case. And 
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Fig. 9. Permeability as function of Dy. 
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- - -- hydrophobicity 


Relative permeability of water 
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Fig. 10. Water relative permeabilities as functions of saturation for different 
hydrophilic pore fraction, d, = 1.2507 (base case). 


water relative permeability of the hydrophobicity case increases 
quickly from 0 to nearly 1 with saturation changing from 0 to 
0.4. But for the same range, the water relative permeability of 
the hydrophilicity case is really very low until the saturation 
comes to 0.6. It begins to increase quickly. The results also show 
that water relative permeability increases with the reducing of 
hydrophilic pore fraction. 

Liquid water is produced in the cathode catalyst layer by 
the electrochemical reaction of PEMFCs. It must be efficiently 
removed otherwise it will cause “flooding” of the electrode 
which hinders the transport of the reactants. Liquid water 
removal from the cathode catalyst layer can be achieved in two 
ways: flow through the cathode GDL to the cathode gas chan- 
nel; diffusion and convection through the membrane to the anode 
side. For the former case, according to the results shown above, 
hydrophobic carbon paper is preferred for its high water rela- 
tive permeability. As for the case of diffusion and convection 
through the membrane, it will be better if the pressure of liquid 
water in the cathode catalyst layer side is high, due to the initial 
pressure needed for the hydrophobic pores. It also can be seen 
that hydrophobic carbon paper is preferred. The conclusion can 
be drawn that a hydrophobic GDL is better than a hydrophilic 
GDL in the view of liquid water removal on cathode side of 
PEMFCs. 

For anode side, liquid water accumulation in the anode side 
is determined by the transport of water in the membrane, which 
is mainly via diffusion, convection and electro-osmotic drag. 
Two situations should be considered: one is the transport direc- 
tion of water from anode to cathode, the other is the trans- 
port direction of water from cathode to anode. If liquid water 
transports from the anode to the cathode, it may cause water- 
depletion on the anode side. So it is necessary to maintain a 
certain quantity of water on the anode side. From this point 
of view, a hydrophilic GDL is suitable for its low water rel- 
ative permeability. However, if liquid water transports from 
the cathode to the anode, there may be an accumulation of 
liquid water on the anode side, so a hydrophobic GDL is pre- 
ferred. 


1.0 


0.8 


06 $ —— hydrophilicity 
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Fig. 11. Water relative permeabilities as functions of saturation for different 
hydrophilic pore fraction, d, = 1. 
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Fig. 12. Water relative permeabilities as functions of saturation for different 
hydrophilic pore fraction, d = 1.5. 


Figs. 11 and 12 show the water relative permeabilities for 
different tortuosity dimensions. The results show that the water 
relative permeability of the hydrophilicity case decreases with 
increase of the tortuosity dimension, and the water relative per- 
meability of the hydrophobicity case increases with increase of 
tortuosity dimension. 


4. Conclusion 


In this study, a fractal model is developed to predict the per- 
meability and liquid water relative permeability of the GDL 
(TGP-H-120 carbon) in PEMFC s for the hydrophilicity case 
and the hydrophobicity case. Area dimension, tortuosity dimen- 
sion, maximum pore size, thickness of GDL, and saturation as 
well as the hydrophilic pore fraction in the GDL are involved 
in this model. The area dimension and tortuosity dimensions 
are obtained by the box-counting method based on the micro- 
graphs (in-plane and through-plane) of TGP-H-120. The model 


852 
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is validated by the comparison between predicted results and 
experimental data. The results indicate: 


. The permeability of the GDL increases with the decrease of 


the tortuosity dimension or the increase of the area dimen- 
sion. 

The water relative permeability of the hydrophobicity case is 
much higher than that of the hydrophilicity case. 


. To effectively remove the water on the cathode side of PEM- 


FCs, hydrophobic carbon paper is better than hydrophilic car- 
bon paper. As for the anode side, if water flows from the cath- 
ode to the anode, a hydrophobic carbon paper is preferred, 
otherwise, hydrophilic carbon paper should be applied. 
Water relative permeability increases with increase of 
the tortuosity dimension for the hydrophobicity case but 
decreases for the hydrophilicity case. 
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